Probing the surface states in Bi 2 Se 3 using the Shubnikov-de Haas effect 
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Shubnikov-de Haas (SdH) oscillations are observed in Bi2Se3 flakes with high carrier concentration 
and low bulk mobility. These oscillations probe the protected surface states and enable us to extract 
their carrier concentration, effective mass and Dingle temperature. The Fermi momentum obtained 
is in agreement with angle resolved photoemission spectroscopy (ARPES) measurements performed 
on crystals from the same batch. We study the behavior of the Berry phase as a function of magnetic 
field. The standard theoretical considerations fail to explain the observed behavior. 

PACS numbers: 73.25. +i, 71.18.+y, 72. 20. My, 79.60.-i 



The recently discovered topological insulators are a 
matter of intense theoretical and experimental study in 
contemporary condensed matter physics (for a review of 
the field, see ref. [l[). The exotic surface states that occur 
in these materials were confirmed and thoroughly stud- 
ied using surface sensitive techniques, such as ARPES 
and scanning tunneling microscopy (STM) Q. These 
measurements demonstrated the absence of backscatter- 
ing at the surface, which is believed to be protected by 
time reversal symmetry. 

Of the 3D topological insulators, Bi 2 Se 3 has the sim- 
plest band structure - a single surface Dirac cone, how- 
ever transport experiments on this material present a 
challenge. The surface conducting state in Bi 2 Se3 is diffi- 
cult to separate from the spurious bulk conductivity, be- 
lieved to originate from Se vacancies 0-0] • Even in the 
lowest carrier concentration samples metallic bulk behav- 
ior is observed. A few successful attempts to surmount 
this problem were reported, including: thinningdown the 
3D material to decrease the bulk contribution [7H10jj , us- 
ing gate voltage 
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pulsed high magnetic fields [13| , 
and changing the material composition to reduce the bulk 
conductivity [El. [lil Il5j|. However, providing conclusive 
transport evidence of the tt Berry phase expected from 
the topology of the surface remains a challenge. 

Recent attempts to probe the surface of Bi2Se3 crystals 
using the SdH effect resulted in an overwhelming bulk 

signal @, EE E3]- These 

groups focused on making the 
samples cleaner in order to increase the bulk resistance. 
However, the decrease in the bulk carrier concentration 
led to an increased mobility. Consequently, the 3D oscil- 
lations dominated the signal. 

In our work we adopt the opposite approach; the 
SdH effect has an exponential sensitivity on the mobil- 
ity, whereas the conductance varies only linearly with it. 
This allows a mobility window in which the SdH effect 
probes only the protected surface states, albeit the dom- 
inance of the bulk over the classical conductance. This 
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FIG. 1: (Color online) Resistance as a function of field for a 
high carrier concentration (SI) at T = 4.2K in a perpendicu- 
lar magnetic field up to 32. 5T. Left inset: typical longitudinal 
resistance versus temperature at H = 0. Right inset: Hall re- 
sistance as a function of magnetic field at T — 2K, H || C3 
axis. The dashed line is a linear fit yielding a charge carrier 



concentration of 5.6 x 10 cm 



can be done by increasing the carrier concentration (thus 
lowering the bulk mobility), in contrast to the conven- 
tional approach. 

Here we report on transport and ARPES measure- 
ments in Bi2Se3 flakes. We observe SdH oscillations in 
the magnetoresistance, and provide evidence that they 
originate solely from the surface Dirac cone. From these 
measurements we find the effective mass and carrier con- 
centration and provide a lower bound for the Dingle scat- 
tering time. The ARPES data provides evidence of a 
single Dirac cone and the Fermi wave vector of the sur- 
face states obtained is in agreement with the magneto 
transport result. We study the evolution of the SdH 
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phase with magnetic field. The broad field range stud- 
ied and the pronounced oscillations observed result in 
unprecedented accuracy at which the frequency can be 
determined. However, we find that the phase is a strong 
function of the frequency. The data is consistent with 
a 7r Berry phase changing with magnetic field. However, 
the conventional Zccman effect alone is insufficient to de- 
scribe our data. Other possibilities, such as zero phase, 
are discussed. 

Single crystals of Bi2Se3 were grown in the standard 
method [i , 18 1 . It is known that a Se rich mixture usually 
results in a reduced carrier concentration in the bulk, as is 
the case in our reference sample with a carrier concentra- 
tion of ~ 10 17 cm~ 3 . In order to get a rather high carrier 
concentration and low bulk mobility, we used a stoichio- 
metric mixture of Bi (99.999%) and Se (99.99%), which 
was cooled down slowly at 3° Cj hour and annealed for 
24h at 620°C. This resulted in a carrier concentration 
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X-ray diffraction spectra (not shown) 



are consistent with previous reports of single crystals of 
Bi 2 Se 3 0. 

Flakes ~ 10— 100/xm thick were freshly cleaved perpen- 
dicular to the C3 axis from a single crystal in a nitrogen 
environment (samples S1-S3). Gold contact wires were 
attached to the sample using silver paint. Measurements 
up to 14T were performed in a Quantum Design PPMS 
platform using Keithley 6221 current-source coupled to a 
Keithley nanovoltmeter 2182A. High field measurements 
(up to 32.5T) were performed in HFML with an SR830 
DSP Lock-In Amplifier. For SdH measurements flakes 
were cut into narrow bars and measured in a four contact 
configuration. Other pieces from the same batch were 
measured in a Van der Pauw configuration for a more 
accurate determination of the resistivity and Hall coeffi- 
cient. ARPES measurements were performed using the 
HeI Q -line (bv=21.2eV) from a Scienta UV-lamp and with 
an R4000 Scienta analyzer. The crystals for the ARPES 
measurement were cleaved in-situ at vacuum better than 
5xl0 _u torr at 20K. 

We took extra care to properly align our voltage leads. 
Misaligned contacts resulted in a spurious Hall contri- 
bution in the longitudinal voltage channel. This contri- 
bution manifested itself in an antisymmetric background 
and a phase shift between positive and negative mag- 
netic fields. Although all samples gave similar frequen- 
cies, samples with such Hall contribution were not used 
for the phase analysis. 

In the insets of Figure [1] we show typical transport 
properties of a sample with high carrier concentration. 
The resistivity as a function of temperature for a typical 
flake is shown in the left inset. A metallic temperature 
dependance with a saturation value of ~ 280/xf2cm be- 
low 20K is observed. The carrier concentration of this 
flake inferred from the Hall resistance at 2K (see right 
inset) is 5.6 x 10 19 cm -3 , this gives a bulk mobility of 
400cm 2 /V ■ sec. Other samples in the literature have 
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FIG. 2: (Color online) (a) Fourier transform (FFT) of the 
oscillations in Figure [1] after substracting a smooth polyno- 
mial background (solid blue) . The dashed red is FFT of the 
same oscillations after dividing out RtRd - see Eq. [3l [4j with 
m* = 0.16m e and to = (3.6 ± 0.4) x 10 _14 sec found in the 
text, (b) Frequency of the oscillations measured as a function 
of 9 (green diamonds - SI, red circles - S2). The bars repre- 
sent the Full Width at Half Maximum of the FFTs, which are 
an upper bound for the error. The solid line is the expected 
2D —s behavior. 



lower carrier density and higher mobilities 1,0. These 
studies also show samples with higher carrier concen- 
tration r~j 10 19 cm -3 , however, their bulk mobility is at 
least a factor of two higher. Apparently, all these re- 
ported mobilities are high enough to obscure surface os- 
cillations. Moreover, our reference (non stoichiometric) 
samples with n ~ 10 17 cm -3 have similar properties as 
other samples in the literature [l6| exhibiting 3D oscilla- 
tions. 

The resistance versus magnetic field at 4.2K is shown 
in Fig. Q] Strong oscillations are observed. Below we 
analyze the oscillations and show evidence that they arise 
from the topological surface states. 

According to the Onsager relation , 
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the frequency of the SdH oscillations F is proportional 
to A(€f), the cross section of the Fermi surface in the 
plane perpendicular to the applied magnetic field, h is 
Planck's constant and e is the electron charge. For the 
case of a two-dimensional (2D) Fermi surface, F oc —^-5, 
where 9 is the angle between the normal to the 2D plane 
and the direction of the magnetic field. In Figure [2] we 
show that our data follows the expected behavior of a 
2D Fermi surface. Here 9 is the angle between the C3 
axis and the applied magnetic field. For 9 = 0, we find 
F = 198, 190, 199T for S1-S3 respectively. 
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FIG. 3: ARPES band dispersion of a high carrier concentra- 
tion sample. Both the bulk and the surface states are observed 
in the data. The dashed lines mark the bottom of the 3D band 
(around 160meV below the Fermi-level) and the Dirac-point 
(around 345meV below the Fermi- level). The arrows mark 
k F of the 2D surface state (k F = 0.084 ± 0.005A" 1 ). 
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FIG. 4: (Color online) Temperature dependence of the SdH 
oscillations at 9 — after subtracting a smooth polynomial 
background (S3). Top inset shows the temperature depen- 
dence of the amplitude, yielding m*~0.16m e . Bottom inset: 
a Dingle plot resulting in td = (3.6 ± 0.4) x 10 _14 sec. 



The behavior shown in Figure [2] is in strong contrast 
with measurements of our low carrier concentration sam- 
ples. For these flakes, F = 20,26T for H || C 3 and 
H _L C3, respectively. This corresponds to a 3D car- 
rier density of 6.6 x 10 17 cm -3 [2(|, consistent with Hall 
measurements. 

Substituting the average frequency F = 196T in Eq. 
[1] and assuming a circular Fermi surface, we find kp = 
0.077±0.003A , which corresponds to a 2D carrier den- 

sity of n2D — 4^ — 4.7 x 10 12 cto~ 2 , assuming no spin 
degeneracy, as should be the case for a topologically pro- 
tected surface. 

Fig. [3] presents ARPES data from a cut taken along 
the TM direction and going through the V point on a 
typical flake. The contribution of the bulk-electrons as 
well as the Dirac-like surface state are clearly seen. We 
find kp — 0.084 ± 0.005A in agreement with the value 
extracted from the SdH analysis. 

The resistance of a 2D system exhibiting SdH oscilla- 
tions is given byfli-j 
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R S 2D = R a {l + R t R d cos[2-k{ 
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where Rq is the zero field resistance and 2ir/3 is the Berry 
phase which is expected to be 7r for an electron rotating 
around the Dirac point in a topological insulator at the 
low magnetic field limit lj . Rp contains the temperature 
dependance, 
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The dingle factor is 



Rd = exp( ) 



(4) 
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where Tp> is the dingle scattering time, and corresponds 
to the dephasing of the Landau states. uj c = is the 
cyclotron frequency. Rp determines the amplitude decay 
with the decrease of magnetic field. 

The standard analysis of SdH oscillations yields a cy- 
clotron mass of to* ~ 0.16to 6 [23| (see Fig. 2]), this to- 
gether with the field dependance of the oscillation ampli- 
tude at 4.2K gives t d = (3.6 ± 0.4) x 10" 14 sec, a Dingle 
temperature of Td = g^nhr = 33.78 ±0.33X and a cor- 
responding quantum mobility of ^ 2 q D = 400cm 2 /V ■ sec. 

Further confirmation of the 2D nature of the oscilla- 
tions reported here is inferred from a comparison between 
the quantum (Dingle) and transport mobilities. Forcing 
a 3D fit to the amplitude of the oscillations (which con- 
tains an additional \f~B coefficient), one finds a 3D Dingle 
time of 3.9 x 10~ 14 sec which corresponds to a bulk quan- 
tum mobility of ^ q D = 570cm 2 /V ■ sec, larger than the 
transport mobility /i^P = 400cto 2 /V ■ sec. This is in 
strong contrast with recent reports for Bi2Se3, where the 
Dingle mobility is four times smaller than the transport 
one [l7l|. Furthermore, it is physically impossible to have 
H q D > fif^ in a bulk material, since the Dingle mobility 
takes into account all scattering events including small 
angle scattering that usually do not affect the resistivity, 
whereas for the transport mobility, backscattering events 
play the major role. 

One of the hallmarks of a topological insulator is the 
7r Berry phase associated with the Dirac dispersion re- 
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lation. Since pronounced oscillations are observed in a 
broad field range (see Figure [lj , we are able to follow 
the evolution of the Berry phase as a function of applied 
magnetic field for SI. 

The standard phase analysis is done by plotting F/B 
versus n, the oscillation index (usually referred to as the 
Landau level fan diagram). Applying this to our data, 
using the FFT peak F — 198T yields a zero phase (ft = 
0). We note that due to the many oscillations measured, 
the Fourier transform analysis favors a frequency that 
yields a zero phase. Furthermore, we find that ft is very 
sensitive to F, and even a deviation of 1.5%, which is well 
within the error margin, has a strong effect on ft. For 
example, if we take F — 195T, we get a ir Berry phase 
(ft = 0.5) at low fields, deviating towards zero as the field 
increases. In the following, we study the behavior of ft for 
frequencies in the vicinity of the FFT peak, F = 198±3T. 

In order to focus on the field and frequency dependen- 
cies of the phase, we use the relation 



ft(F,n) =n- F/B 



(5) 



and plot ft versus n in Figure [5l where n is the index 
of the oscillation minima for the two frequencies men- 
tioned (see supplementary information for more details). 
The two behaviors mentioned above are now clearly seen. 
This analysis casts strong doubts on the ability to inde- 
pendently determine the Berry phase using SdH oscilla- 
tions. 

Assuming that at low field ft should be 0.5, the sim- 
plest explanation for the phase change would be a Zee- 
man term in the Hamiltonian [l3| . In Figure [5] we plot 
the expected behavior of ft using gyromagnetic ratios 
g = 50, 100 and vp = 6.24 x 10 5 m/sec as inferred from 
the ARPES measurement (see Figure [3]). While both 
theoretical curve (g — 100) and the measured data are 
in qualitative agreement, it is obvious that this simple 
model is insufficient to describe our measurements. It 
is possible that the non-ideality of the Dirac dispersion 
should also be taken into account 21 j . 

In summary, we present evidence that the surface 
states in Bi 2 Se 3 can be probed in highly conducting 
flakes using the Shubnikov-de Haas (SdH) effect. This 
is in contrast with the conventional approach focused 
on improving the crystal purity. From the SdH analy- 
sis we find for the surface states m* = 0.16m e , tjj — 
(3.6 ±0.4) x 10 _14 sec. The Fermi momentum is in agree- 
ment with the ARPES data obtained on a flake from the 
same batch. We carefully study the behavior of the Berry 
phase with magnetic field, and show that two scenarios 
are possible: ft = 0.5 (ir Berry phase) at low magnetic 
fields which changes fairly quickly with magnetic field, 
or a trivial zero phase. Two major issues are yet to be 
understood: first, the survival of protected surface states 
despite the large bulk conductivity and second, the pe- 
culiar behavior of the Berry phase with magnetic field. 
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FIG. 5: (Color online) Analysis of the Berry phase 2nft ac- 
cording to Eq. [5] using the data from Figure [1] after subtract- 
ing a smooth polynomial background, n = 6.5 corresponds 
to B = 30.4T and n = 28 to B = 7T. The blue diamonds are 
obtained using the FFT peak of F = 198T, consistent with a 
zero Berry phase. The black circles correspond to F = 195T, 
for which a tv Berry phase is observed at low fields. This phase 
changes towards zero as the field increases. This demonstrates 
the strong sensitivity of the phase to the frequency chosen. 
The error bars stem from the freedom in choosing the oscilla- 
tion extrema. The green dashed (red dashed-dot) line is the 
expected behavior of ft taking into account a Zeeman term 
with g — 50 (100) respectively. 
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